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ABSTRACT 

Sulfur appears to be depleted by an order of magnitude or more from its elemental 
abundance in star-forming regions. In the last few years, numerous observations and 
experiments have been performed in order to to understand the reasons behind this 
depletion without providing a satisfactory explanation of the sulfur chemistry towards 
high-mass star-forming cores. Several sulfur-bearing molecules have been observed in 
these regions, and yet none are abundant enough to make up the gas-phase dehcit. 
Where, then, does this hidden sulfur reside? This paper represents a step forward in 
our understanding of the interactions among the various S-bearing species. We have 
incorporated recent experimental and theoretical data into a chemical model of a hot 
molecular core in order to see whether they give any indication of the identity of the 
sulfur sink in these dense regions. Despite our model producing reasonable agreement 
with both solid-phase and gas-phase abundances of many sulfur-bearing species, we 
find that the sulfur residue detected in recent experiments takes up only ^6 per cent 
of the available sulfur in our simulations, rather than dominating the sulfur budget. 

Key words: astrochemistry - ISM: abundances - ISM: molecules - solid state: 
refractory - solid state: volatile - stars: formation. 


1 INTRODUCTION 


Sulfur, despite being only the tenth most abundant element 
in the Milky Way, is of significant astrochemical interest. 
Two of the key open questions regarding sulfur are: (i) where 
is the sulfur that appears to be depleted from the gas phase 
in dense regions? and, (ii) what use can studies of sulfur¬ 
bearing species be in our understanding of astronomical en¬ 
vironments? 

The depletion of sulfur became evident in the 1970s, 
80s and 90s, when the observed abundances of sulfur¬ 
bearing molecules in dense regions did n ot match the 
observed cosmic abundance of sulfur (e.g., Pemiaset_^lJ 


197ll: ODPenheimer fc Dalgarnd[l97l : lTieftrunk et al.lll994 


Palumbo et al.lll997l l. whereas in diffuse and highly ionised 
regions, abu ndances of sulfur seemed roughly cosmic 
(-^10~^; e.g., Savage fc Sembac^ 19961: Howk et all 200d; 


iMartm-Hernandez et al.l l2002l : iGarci'a-Roias et al.l 1200^ . 
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This is often referred to as the ‘sulfur depletion problem’. 
Since then, the interest in the chemistry of sulfur has height¬ 
ened. To comprehend the sulfur depletion problem fully it 
is important to go through the results that have been pub¬ 
lished on this topic, in order to appreciate the state-of-the- 
art; especially since relevant experimental work has been 
performed recently. 

In 1999, it was proposed that sulfur typically existed 
in ionised form (S^) in translucent gas, and thus froze 
out more rapidly than neutrals upon cloud collapse, due 
to electrostatic att raction to the negatively-charged grains 
llRuflie et al.lll999ll . However, the form which sulfur takes 
upon freeze-out is not evident: gas-phase carbonyl sulfide, 
PCS, has been detect ed in star-forming regions (SFRs; 
Ivan der Tak et ^l2003l l. at a level one thousand times too 
low to be the dominant carrier of sulfur; but, estimates of 
grain-surface PCS abundances may be affected by blending 
with a nearby methanol feature in infrared spectra. Chem¬ 
ical models at the time suggested that S, SP, CS and H 2 S 
may be viable sinks of sulfur atoms (e.g., iMillar fc HerbstI 
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I 1990 I : Ijansen et alJ Il995l ~) , but our understanding of sulfur 
chemistry, particularly the chemistry of CS, is not complete. 

Our second key question asks about the intelligence that 
sulfur-bearing species can give us in the understanding of 
astronomical environments. It has been suggested by sev¬ 
eral authors that sulfur-bearing species may act as evolu¬ 
tionary tracers for a specific r egion . In an attempt to study 
this potential role. ICharnlevI (Il997l l proposed that SO/H 2 S 
and SO/SO 2 ratios act as molecular clocks for grain man¬ 
tle disruption since these ratios seem to vary between dif¬ 
ferent astronomical environments (dark clouds, hot cores, 
shocks and winds aroun d protostars were stu died) and also 
within individual SFRs. iHatchell et al.l lll998ll followed this 
idea, to constrain the ages of cores, by looking at the ratios 
of sulfur-bearing species towards ultra-compact HII regions. 
In particular, they developed different chemical models for 
each of the eight sources t hey investigated , by v arying dif¬ 
ferent physical parameters. IHatchell et al.l (Il998l ') classified 
the sulfur-bearing species along with the age of the young 
stellar object: 

• H 2 S, SO are abundant in younger cores 

• H 2 S, SO, SO 2 at intermediate ages 

• later, SO and SO 2 are present, but without H 2 S 

• hnally, CS, H 2 CS and PCS become the most abun¬ 
dant sulfur-bearing species. IHatchell et al.l suggested OCS 
is formed on the grain surface. 

A few years later, IViti et ^ HooJ) proposed that NS/CS 
and SO/CS ratios were specihc indicators of a shock pas¬ 
sage in the vicinity of a hot core. In these physical con¬ 
ditions, the sulfur chemistry was found to be connected 
to the HCO/H 2 CO ratio. High values of these ratios indi- 
cated that a shock ha d passed through the m edium. Finally, 
IWajmlam et al.l (l2004h repeated the study bv IHatchell et all 
1 1998ll and highlighted how none of the ratios involving the 
four most abundant sulfur-bearing species (H 2 S, OCS, SO 
and SO 2 ) could be useful by itself for estimating the core 
ages, because the amount of each molecule depends at least 
on the physical conditions, the adopted grain mantle compo- 
siti on and also evolutionary stag e . A re latively recent paper 
by IWakelam. Hersant fc HerpinI |20l3) reported the study 
of S-bearing species in four different high-mass dense core 
sources in order to invest igate the depend ence of their abun¬ 
dances along with time. IWakelam et al.l were unable to re¬ 
produce the observed abundances for OCS, SO, SO 2 , H 2 S 
and CS, but they found that the ratios between OCS/SO 2 
and H 2 S/S O 2 could be used to constrain some evolutionary 
time-scales. I Wakelam et"^ also highlighted the difficulty in 
reproducing the amount of CS, which was overestimated due 
to the fact that its abundance varies with radius and there 
is an uncertainty about the location of the emitting region. 

In order to fully understand what we know to date 
about the presence of sulfur-bearing species in the inter¬ 
stellar medium, we present a summary of all the species 
belonging to this family of molecules which have been ob¬ 
served either in the gas-phase (Table [T]) or on the grain 
surface. In the gas phase, sulfur is a ubiquitous element: 
it has been detected in different astronomical environ- 


ments, from the diffuse medium (|Lisz ^ l2009ll to dark clouds 


llDickens. Laneer & Velusamvl 

200fllL as well as in hot cores 

and hot corinos (ISchoier et al. 

I 2 OO 2 I: ISutton et al.lll995h: in 


comets teoissier et al.ll2007^ : in evolved stars Jwoods et aP 


I 2 OO 3 II ; and in the atmosphere of Venus llKrasnopolskvilioOSl l , 
in various chemical forms. Its emission, therefore, has been 
widely observed, and that has enabled the depletion of sul- 
fur to b e stud ied in a variety of environments. For instance, 
IjenkinsI (l2009lj observed atomic sulfur lines towards the dif¬ 
fuse medium and showed a relationship between the amount 
of depletion of the elements and the density of the cloud. 


On the other hand, the only two S-bearing species 
Hrmly detected on grain surfaces have been OCS, 
with relatively l ow fractional abundances o n the 
order of 10~^ dPalumbo. Tielens fc Tokunagal 19951: 
Palurnbo. Geballe fc Ti elensI 199lt l and SOt i Boogert et alJ 

and thus the form of grain- 


1993; IZasowski et al.l 


surface sulfur is relatively unknown. Several experiments 
were performed in order to understand which molecules are 
candidates for e xplaining sulfu r chemi stry in the solid state. 
In pa rticular, iFerrante et al.l l|2008h and ICarozzo et al.l 
(|20inh investigated the mechanism of formation of OCS, 
discovering that CO reacts with free S atoms produced by 
the fragmentation of the sulfur parent species. OCS was 
seen to be readily formed by cosmic-ray irradiation, but 
at the sa me time it was easily destroyed after continued 
exposure. IFerrante et al.l l|2008l ~) observed CS 2 production 
as one of the main product channels, although carbon 
disulfide has not been yet detected in interstellar ices. An 
alt ernative molecule, hy drated sulfuric acid, was suggested 
bvIScap pini et a P ll200,‘^ as the m ain sulfur reservoir. Later 
iMoore. Hudson fc CarlsCT |2003) produced this species 
by ion irradiation of SO 2 and H 2 S in water-rich ice over 
the temperatu re range 8 6-130 K. More recent papers by 


ICnrozzo et ^ (l201Clli; Jime nez-Escobar fc Munoz Card 
( 2011^ and I.Iimenez-Escobar. Munoz Caro fc ChenI ( 2014h 
focused on the products of cosmic-ray and UV-photon 
irradiation of H 2 S ice analogues. Among the products, they 
found the presence of a sulfur residue which might explain 
the missing sulfur in dense clouds. A revised chemistry of 
S-bea ring species has been discussed bv IPruard fc WakelamI 
(|2012tl . who were able to explain, in an initial attempt, the 
lack of observation of H 2 S on the grain surface, but without 
considering the presence of refractory sulfur. 


The above summary of the state-of-art shows that, de¬ 
spite several investigations, sulfur chemistry in the inter¬ 
stellar medium seems an intricate problem that is yet to 
be solved. The present paper is a step forward in our under¬ 
standing of the sulfur chemistry in regions of star-formation. 
In particular, we couple recent experimental data with a re¬ 
vised version of the ucl_CHEM chemical model, taking into 
account the formation of this sulfur residue. Moreover, we 
present a new classification of the thermal desorp tion of sul¬ 
fur sp ecies compared to the one elaborated by IViti et ^ 
ll2004h . Non-thermal desorption effects are also taken into 
account. The paper is organised as follows: Section [5] de¬ 
scribes the astrochemical model; Sect. |3] contains the de¬ 
tails of all the experimental and theoretical results that 
have been included in our chemical modelling, described 
in two different subsections; Sect. St] contains the outputs 
from UCL_CHEM models with their analyses and discussion; 
Sect. 14.21 compares our theoretical results with data from 
the observations; and finally, we present our conclusions in 
Sect. [ 5 ] 
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Table 1. List of gas-phase sulfur-bearing species with the dense sources towards which they were first observed. Example references are 
provided in the third column. 


Molecule 

Source 

Column density [cm 

cs 

Orion A, W51 

2 - 210 x 10 ^^ 

DCS 

SgrB2 

^SxlQi® 

H 2 S 

Hot cores 

4-50x101® 

H 2 CS 

SgrB2 

>lxl 0 i® 

SO 

Orion A 

~lxl0i5 

SO 2 

Orion A 

3-35x10®® 

SiS 

IRC-l-10216, SgrB2 

4x10®® 

NS 

SgrB2 

1 x 10 ®® 

CH 3 SH 

SgrB2 

2 x 10 ®® 

HNCS 

SgrB2 

3x10®® 

HCS+ 

SgrB2, Orion 

2 - 200 x 10 ®® 

C 2 S 

TMC-1, SgrB2, IRC-l-10216 

6-15x10®® 

C 3 S 

TMC-1, SgrB2 

1 x 10 ®® 

SO+ 

IC 4434 

~5x10®2 

HSCN 

SgrB2(N) 

1 x 10 ®® 

SH+ 

SgrB2 

< 2 x 10 ®® 

HS 

W49N 

5x10®® 

CH 3 CH 2 SH 

Orion KL 

2 x 10 ®® 


Reference 

^enzia^e^ah 

Jeffert^^al^ 

^haddeu^e^aL ^^72) 

^inclai^^L ^^73 ) 

Gottlie^^^B^ ) 

^n^de^^aL ^^75) 

^Morris^^t_^ (1975) 

Gottlie^^^L ; JOu^e^^aL ^^7^ 

^^nkej_^terkin^^^_Thaddei^ (1979) 
FYerking^Link^^^Thadde^ ^^79) 
^haddeus^Guelii^^^inl* ^^^) 

^^8^ ; ^ernichax^^^L ^^8^ 
I^aib^^aL ^^^) ; Yamamot^^^ah ^^^) 
^^irner_et_^ (1992) 

Halfei^^aL ^009) 

Mentei^^L ) 

Neufel^^^L 
Kolesnikova et al. (20141 


2 MODELLING 


We m odify a pre-existinR model, the ucl_C HEM chemical 
code dViti fc William j Il999l : IViti et aD l2004l l , in order to 
include the new experimental results. Before going into the 
detail of our updates, we briefly describe the physics and 
the chemistry behind the model. The model performs a 
two-step simulation. Phase I starts from a fairly diffuse 
medium where most chemical species are in atomic form 
(apart from H 2 ), which undergoes a free-fall collapse un¬ 
til densities typical of the gas that will form hot cores or 
hot corinos are reached (10^ cm“® and 10® cm“®, respec¬ 
tively). During this time, atoms and molecules are depleted 
on to the grain surfaces and they hydrogenate when possi¬ 
ble. The depletion efficiency is determined by the fraction 
of the gas phase material that is frozen on to the grains. 
This approach allows a derivation of the ice composition 
by a time-dependent computation of the chemical evolu¬ 
tion of the gas-dust interaction process. The initial elemen¬ 
tal abundances of the main species (such as H, He, C, O, 
N, S and Mg) are the main inputs f or the chemistry, and 
are taken f rom [Sof ia &: M ever booill and other references 
detailed in IViti fc Williams! (Il999ll . in common with much 
recent ucl_CHEM work. We usually assume that, at the be¬ 
ginning, only carbon and sulfur are ionised and half of the 
hydrogen is in its molecular form. The other elements are 
all neutral and atomic. Gas-phas e reactions are taken from 
the UMIST RATE06 database J Woodall et al.l [20011 1. and 


freeze-out reactions are included in the reaction network in 
order to allow the formation of mantle species. The deple¬ 
tion efficiency can be modified by adjusting the freeze-out 
fraction parameter fr. In Phase I we follow both the gas- 
phase and the grain-surface chemistry, and the transitions 
between these two phases (freeze-out and desorption). 


Phase II is the warm-up phase and the model follows 
the gas-phase chemical evolution of the remnant core when 
the protostar itself is formed. During this stage, the impor¬ 
tance of grain-surface reactions decreases with the increas¬ 
ing temperature, because of the sublimat ion of importan t 
molecules (such as CO) even at ~20K (see lViti et al.ll2004h . 


The time-dependent evaporation of the ice is treated in one 
of two ways. Before running the phase II model, we can 
choose a final gas temperature (maxt) for the astronomi¬ 
cal object studied. The treatment of evaporation can be ei¬ 
ther time-dependent (where mantle species desorb in various 
temperature bands according to the experimental results of 
ICollings et al.lf2004l or instantaneous (in that all species will 
desorb from the grain surfaces at the hrst timestep). In this 
paper we only consider time-dependent thermal desorption. 
Non-thermal desorption of spe cies is also taken into account 
and is based on the study by [Roberts et all ll200it) . Three 
desorption mechanisms are included in ucl_chem: desorp¬ 
tion resulting from H 2 formation on grains, direct cosmic-ray 
heating of the ice, and cosmic-ray photodesorption. The lat¬ 
ter mechanism is due to the generation of UV photons which 
occur when cosmic-ray particles impact the grain. We do not 
account for direct UV photodesorption, since the density, 
and hence UV extinction, in these regions is generally high. 

Finally, the outputs from the code consist of the frac¬ 
tional abundances of all the species (in both gas phase and 
on the grain surface) as a function of time. Fractional abun¬ 
dances of molecules are calculated as a ratio to the total 
number of hydrogen nuclei (n(H) -|- 2n(H2)), where n rep¬ 
resents the number density in cm“®. 


3 UPDATE OF UCL.CHEM WITH RECENT 
RESULTS 

3.1 Experimental results 

We have updated the ucl_CHEM gas-grain chemical model 
to include some published and unpublished experimental 
data on reactions involving S-bearing species occurring on 
icy mantle analogues. Our aim is to derive more informa¬ 
tion about the existence of a sulfur reservoir. First of all, 
we summarise the results obtained from the different exper¬ 
iments (see Subsections [ 3 m [3T^ that we insert into the 
chemical model. Then, we run models to benchmark the dif- 
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Figure 1. Column density of selected species as a function of ion 
fiuence after irradia tion of a CO:H 2 S = 10:1 ice mixture at 20K, 
using the data from iGarozzo et alj ll201Cl i. Solid and dotted lines 
have been drawn to guide the eye. 

ferences that the new data make on the sulfur chemistry (see 
Sect. 2]). 

3.1.1 The effect of cosmic rays on icy H 2 S mantle 
analogues 

Cosmic rays impinging on icy mantles are energetic enough 
to induce chemical and structural modifications on the grain 
surface. As a consequence the ice constituents differ from 
the composition of the gas, and experiments of astronomical 
relevance are therefore the only tools that can provide us 
with information concerning the potential interactions that 
can arise after these dynamic impacts. 

In order t o sim ulate a flux of cosmic-ray particles, 
iGarozzo et al.l ll2010ll irradiated their ice sample with 200 
keV protons at 20 K in a high-vacuum chamber (P ^ 
10“^mbar). Their sample consisted of a CO:H2S=10:1 mix¬ 
ture. IR spectra were recorded before and after irradiation in 
order to estimate the actual amount of each species produced 
from the ice mantle analogues due to the proton bombard¬ 
ment. The integrated intensity measured for each selected 
band (in optical depth t{v) units) is proportional to the col¬ 
umn density of the species itself. Column densities for each 
species, N(X), were calculated as the ratio with respect to 
the initial amount of H 2 S, Ni(H 2 S), in the mixture. Figure[T] 
shows the decrease of the initial amount of hydrogen sulfide 
as the products of its dis sociation form. The failure to detect 


H 2 S 

in the solid phase (lEhrenfreund. Charnlev & Wooden 

2004 

: ICarozzo et al.ll2010l: Ijimmez-Escobar & Munoz Caro 

2011 

) may therefore be linked to the strong (~ 80 per cent) 


reduction in its column density when it is sub jected to irra¬ 
diatio n. This is in fact what was postulated bv icodella et al.l 
(l2006li in order to justify the presence of a large amount of 
gas-phase OCS observed in an extended, high velocity gas 
in the massive SFR, Cep A E ast. We have a n alyse d the ex¬ 
perimental data collected by ICarozzo et aP (l2010l l. and in 
particular, we have fitted their data with an exponential 


Figure 2. Column density ratio N(H 2 S)/Ni(H 2 S) of hydrogen 
sulfide as a function of ion fiuence (>1?) after irradiation of th e ice 
mixture. Data points are taken from ICarozzo et al.l ll2010l i and 
fitted with an exponential curve. 


Table 2. A list of the experimentally-detected sulfur-bearing 
species and their laboratory production cross sections. 


Molecule 

p [cm^] 

H 2 S 2 

l.OSxlO-^'^ 

SO 2 

6.78x10-1® 

CS 2 

2.64x10-1® 

OCS 

2.05x10-1® 

S-residue 

1.19x10-1® 


curve (as plotted in Fig. [2} in order to evaluate the reaction 
cross section (a— 4.7x10“^® cm^). The identified molecules 
with their specific production cross sections are listed in Ta¬ 
ble [ 2 I The latter parameter was extrapolated by fitting the 
experimental data at low fiuence with a straight line. 

Three important c onsiderations c o me fr om the analysis 
of the data taken from lCarozzo et al.l (1201011 : 

(i) A detection of CS 2 has not yet been made in the ISM. 
Laboratory data allow us to quantify the abundance of this 
molecule and to extend our model to include a reaction 
scheme for its formation and loss channels. 

(ii) There is evidence of a residue of species containing 
sulfur, that can be estimated as follows: 

iV(S-residue) = Ni{n 2 S) - iV(H 2 S) - A'(S02) - Af(OCS) 

- 2A(CS2) - 2iV(H2S2) (1) 

This residue provides a constant supply of sulfur and it may 
affect the rol e of sulfur as an e volut ionary tracer. Moreover, 
as stated by lAnderson et ^ ll2013ll . the sputtering of this 
residue from the surface could lead to the release of the 
large amount of atomic sulfur seen in shock regions. 

(iii) These experiments help to predict the different forms 
that sulfur can take on the grain surface and to estimate the 
ratios among sulfur-bearing molecules. 
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Table 3. Dissociation of solid H 2 S due to cosm i c-ray impact (CR) 
experimentally investigated bv iGarozzo et aT UmO). All the re¬ 
action channels are provided with a rate (in s“^) of ISM relevance. 
The m before the molecular formula stands for mantle. 


Reaction 



/c = a [s ■*■] 

mH2S + 

CR 

mSres + H 2 

9.53x10-1^ 

mH2S -b 

mCO 

mOCS -b H 2 

1.65x10-1’’ 

mH2S -b 

mH2S mH2S2 -b H 2 

8.23x10-1’’ 

2 mH2S -b 

mCO 

—>■ mCS 2 -b mO -b 2 H 2 

2 .12x10-1^ 

mH2S -b 

2mO 

— y mS02 H 2 

5.42x10-1’’ 


We therefore insert into UCL_CH EM the chemica l reac - 
tions experimentally investigated by iGarozzo et al.l (120101 1 
and reported in Table |3] Note that the reactions listed in 
Table m are simplified representations of a complex exper¬ 
imental process. The action of the high-energy protons on 
the laboratory ice analogues causes ~10® molecular bonds 
to break, leading to a chain of recombination reactions. As 
a result of these very rapid recombinations, the products on 
the right-hand side are observed. These reactions do not take 
place in thermodynamic equilibrium. The remaining sulfur, 
or mSres, stays as a refractory element on the surface. 

The rate constants, k (s“^), have been evaluated by the 
product of the reaction cross section {a in cm^) and the 
flux of cosmic ions (Tism in cm“^s“^). To apply the labo¬ 
ratory results to interstellar medium conditions both quan¬ 
tities have been corrected using the following assumptions: 

(i) We derive the flux of cosmic ions in th e approximation 
of eff ectively mono-energetic 1 MeV protons (iMennella et al.l 
I 2 OO 3 II : Tism = 0.22 ions cm ^s Fism must be regarded 
as an effective quantity: it represents the equivalent flux 
of 1 MeV protons which gives rise to the ionisation rate 
produced by the cosmic ray spectrum if 1 MeV protons 
were the only source for ionisation. 

(ii) Furthermore, we speculate that the cross section 
scales with the stopping power (SP, the energy loss per 
unit pat h length of impinging ions) . According to the SRIM 
code by IZieeler fc BiersacM ( 2009li , in the case of protons 
impinging on a CO:H 2 S mixture, 5P(200keV protons) = 
2.7x5P(lMeV protons). 

In the case of hydrogen sulhde we define (tism = cr/2.7 cm^ 
and we derive a destruction rate of hydrogen sulhde equal to 
3.8x10“^® s“^. In the case of species formed after irradiation 
we dehne ctism = p/2.7 cm^ (see Table [2] for p values) and 
the calculated formation rates for each species are given in 

Table 0 

In our code two-body grain-surface reactions are consid¬ 
ered to be bimolecular reactions occurring as they would in 
the gas phase, meaning that the parameter a should be ex¬ 
pressed in units of volume (cm® s“^). Since the data from ex- 
peri ments are in units of time , we have followed a procedure 
(see lOcchiogrosso et al.ir201ll l for evaluating the rate coefh- 
cients for the formation of S-bearing species in the icy man¬ 
tle, viz-a-viz, we consider an excess of one of the reactants. 
The rate of the reaction will therefore vary only with the 
concentration of the second reactant. As the total amount 
of ice varies with time, we calculate a value for the rate of 
each reaction that varies with time. For instance, when we 
consider CO as the most abundant reactant (this is the case 


of the second and the fourth reaction in Table 0 , reaction 
rates cover six orders of magnitude because of the wide in¬ 
terval spanned by the CO abundances during the collapse 
phase of protostar formation. We point out that density and 
freeze-out rate indeed play a pivotal role in controlling the 
trends of the molecular abundances and since both of these 
two physical parameters only significantly change towards 
the end of the collapse, we therefore refer only to the final 
abundance of CO in order to scale our reaction rates. 

3.1.2 Laboratory investigations of solid OCS formation 

In the experiments at the Cosmic Dust Laboratory (Depart¬ 
ment of Chemistry, University College London) reactants 
are co-deposited on to a highly ordered pyrolytic graphite 
(HOPG) substrate. During the dosing period, the substrate 
is held at a constant temperature, typically in the range 
12-100 K. Once the reactants are deposited, the material is 
allowed to cool to 12 K. Following the deposition, the sample 
is heated up to 200-300 K (depending on the system stud¬ 
ied) during which time mass spectra are recorded. The result 
is a histogram of ion intensity as a function of the ion mass- 
to-charge (jn/z) ratio and the surface temperature. These 
experiments are repeated at different temperatures to give 
the dependence of the amount of product formed on the sub¬ 
strate during dosing. Finally, a simple kinetic model is run to 
derive reaction barriers and Arrhenius pre-exponential fac¬ 
tors from the temperature pro file, allowing rate constant s 
to be calculated. In particular, IWard. Hogg fc Pried (l2012li 
studied the following route for the formation of solid OCS: 

mCS 2 -f mO ^ mOCS -b mS, (2) 

with CS 2 and O that are both adsorbed on 
the icy mantles. They found a rate constant of 
1.24x10“^® cm^ molecule”^ s“^. As mentioned in the 
previous section, since our code accounts for reactions that 
occur in three dimensions (i.e. as gas phase reactions), 
we need to tran sform the surface rate constants from 
IWard et al.l (l2012l l into standard units, cm® s“^ (for more 
details o n the theoretical assump tions in this transforma¬ 
tion, see lOcchiogrosso et al.ll2012ll . We therefore calculate 
a final value of the rate constant as 6.4x 10“^® cm® s“^ at 
20 K for the reaction between O and CS 2 . 

3.2 Incorporation of new theoretical data into the 
model 

In order to investigate the form of sulfur once it freezes on to 
the grain surface, we insert an extended chemistry including 
all the S-bearing species mentioned in the previous Sub¬ 
section into our gas-grain chemical network. In addition to 
these reactions, we also insert two new paths for the forma- 
tio n of carbonyl s u lfide, OCS, as theoretically investigated 
by I Adriaens et ahl (l201Cll L Reactions are listed in Table |4] 
The rate parameters relative to each channel are also 
reported; 7 represents the reaction barrier in units of Kelvin 
(K). Note that, unlike lAdriaens et al.l (l2010lL we do not ac¬ 
count for the formation of any adduct species and we do 
not distinguish among the cis-trans geometries of the reac- 
ta nts and the products. The adsorption energies calculated 
by lAdriaens et al.l 1 I 2 OIOI ') were smaller in value than those 
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Table 4. Rout e s to P CS formation on a coronene surface, from 
lAdriaens et alJ l|2QlQh . a, /3, 7 represent the parameters for the 
rate coefficient in the modified Arrhenius equation. They have 
been adapted for a water-ice surface by Adriaens (2013, priv. 
comm.). 


Reaction 

Q [S -I] 

/3 

7 [K] 

CO + S ^ OCS 

1 .66x10“^^ 

0 

1893 

CO -1- HS ^ OCS -1- H 

1 .66x10-11 

0 

831 

Table 6. Physical parameters 

for the model 

of a 

prototypical 


high-mass star, fr gives an indication of the freeze-out efficiency, 
maxt is the maximum temperature reached in the model, size is 
the diameter of the hot core, dens (df) is the (maximum) density 
reached in the collapse, at a time earlier than tfin. represents 
the cosmic ray ionisation rate. 


Parameter 

Value 

temp [K] 

10 

fr [%] 

98 

maxt [K] 

100 

size [pc] 

0.03 

dens [cm“^] 

2 x10® 

df [cm-®] 

IxlOl 

tfin [yr] 

10 " 

C [s-l] 

1.3x10-11’ 


experimentally determined in previous studies l|Piper et al.l 
1 19841 : iMattera. et al.l llQSOl ) due to the fact that the theory 
considers a perfect coronene surface and neglects the weak 
physisorption interactions observable in the presence of sub¬ 
strate defects. Every reaction activation barrier given can 
therefore be seen as an upper limit to the effective barrier. 

In addition to incorporating a grain-surface reaction 
network for sulfur-bearing species as a revision to the orig¬ 
inal version of the code, we also update our gas-phase re¬ 
action network with new rate coefficie nts taken from the 
KID A database dWakelam et al.l I 2 OI 2 I ). In particular, we 
insert or amend the rate parameters from the UMIST 
d atabase for dif f erent OCS routes of formation as described 
in lLoison et al. I ll20ll see Table E]). 


4 PRESENTATION OF RESULTS 
4.1 Trends among the data 

We commence by setting the initial fractional abundance 
of gaseous sulfur ions equal t o 1.4x10“® (as measured by 
ISofia. Cardelli &: Savagg[l994h . We also fix the physical pa¬ 
rameters to values typical of high mass SFRs, listed in Ta¬ 
ble El 

Since we aim to investigate the effect of the new sur¬ 
face reactions on the abundances of the S-bearing species, 
we have to allow molecules to freeze on to the grain. In par¬ 
ticular, we consider two limiting cases where each species 
freezes as itself (FI) or it hydrogenates (F2) and an inter¬ 
mediate case (F3) where we hydrogenate half of the accret¬ 
ing S-bearing species. The grid of the freeze-out pathways 
chosen is reported in Table 0 Note that not all the sulfur¬ 
bearing species included in our model are shown; in partic¬ 
ular, we assume that H 2 S, H 2 CS, OCS, SO 2 , NS, SO, CS 2 
and H 2 S 2 accrete on to the grains without hydrogenating. 


Table 8 . Fractional abundances (with respect to H 2 ) of solid 
sulfur-bearing species obtained as outputs from our code before 
(OLD) and after (NEW) our updates at the end of Phase I. The 
m before the molecular formula stands for mantle. 


Species 

OLD 

NEW 

mCS 

2 .6x10-1® 

2 .6x10-1® 

mOCS 

3.9x10-®® 

6 .6x10-0® 

mHCS 

1 .6x10-®® 

1 .8x10-0® 

mH2S 

1.4x10-®® 

1.4x10-0® 

mH2CS 

5.8x10-®® 

5.7x10-0® 

mS 

trace 

trace 

mS2 

5.8x10-1® 

2.4x10-1® 

mSO 

8 .8x10-1® 

9.2x10-10 

mS02 

5.4x10-1® 

5.7x10-10 

mNS 

9.9x10-1® 

1 .1x10-11 

mHS 

trace 

trace 

mHS2 

3.2x10-1"! 

3.7x10-1® 

mH2S2 

1 .1x10-1® 

2 .2x10-10 

mCS2 

none 

8.7x10-1® 

mS-residue 

none 

5.4x10-0® 


In addition, we need to derive an estimate of the sul¬ 
fur species in the gas-phase, after their thermal evaporation 
from the grain surf ace, for comparis on to evolved hot cores. 
As investigated bv IViti et al.l (|2004| ). species can evaporate 
from the i cy mantle i n different bands of temperatures; in 
particular, IViti et ahl determined that H 2 CS behaves as a 
H 2 0-like molecule that co-desorbs with the H 2 0-ice when it 
starts to sublime from the grain surface (at ~100K). More¬ 
over, the authors classified HCS, OCS, H 2 S, SO 2 as inter¬ 
mediate species since they showed two peaks (due to vol¬ 
cano and co-desorption effects) in their Temperature Pro¬ 
grammed Desorption (TPD) traces. After further discussion 
(W. A. Brown, priv. c omm.), we have revised the above des¬ 
orption classification JViti et al.ll2004l ) as follows: 

(i) HS, H 2 S 2 , OCS, H 2 S, SO 2 , HCS, NS as intermediate 

(ii) H 2 CS, SO as water-like 

(iii) HS 2 as reactive 

(iv) S 2 , CS 2 , S-residue as refractory. 

No sulfur-bearing species are classified as having a desorp¬ 
tion be haviour that is CO-like (the initial desorption cate¬ 
gory of IViti et al.ll20oJ ). 

In order to assess the impact of our updates on the 
chemical network, we ran some test models, comparing our 
original version of the code and those including the new 
sets of reactions. Our findings are reported in Tabled The 
changes in the ice composition are more evident at the end of 
the collapse phase (Phase I), when temperatures are as low 
as 10 K. Note that we only show our results assuming the 
highest degree of hydrogenation (F2) during the freeze-out 
of the molecules and we reserve a more detailed investigation 
for later in the paper. 

Scrutinising Table [SI we immediately observe that the 
most abundant sulfur species in the solid state is H 2 S, which 
is predictable since we are discussing the model with the 
highest degree of hydrogenation. The abundance of OCS is 
enhanced in the new model by a factor of ~2; this is mostly 
due to the rea ction of H 2 S + CO on the grain surface, with 
th e rate from iG arozzo et al.l 1 I 2 OIOI ). The reaction studied 
by I Ward et al.1 (l2012l ) contributes at the level of a few per¬ 
cent. The most interesting result is definitely the amount of 
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Table 5. Gas pha se paths of OCS formation and destruction and their competitive routes. The revised a. values are taken from 
iLoison et al H2oi3). 


Reaction 



New a. values 
[cm^mol“^s“^] 

Previous a values 
[cm^mol~^s~^] 

CH + SO 


OCS -1- H 

1 .1x10"^" 

- 

CH + SO 


SH -1- CO 

9.0x10-11 

- 

0 + HCS 


OCS -1- H 

5.0x10-11 

5.0x10-11 

0 + HCS 


SH -1- CO 

5.0x10-11 

- 

H + HCS 


H 2 -f CS 

1.5x10-10 

4.0x10-0 

S + HCO 


OCS -1- H 

8 .0x10-11 

- 

S + HCO 


SH -1- CO 

4.0x10-11 

3.6x10-10 

OH + CS 


OCS -1- H 

1.7x10-10 

g.ixio-i"! 

OH + CS 


SH -1- CO 

3.0x10-11 

- 

C + OCS 


CO + CS 

l.OxlO -10 

1 .6x10-0 

CH + OCS 

-A- 

CO + CS -1- H 

4.0x10-10 

- 

CN + OCS 

-A- 

CO + NCS 

l.OxlO -10 

- 


Table 7. Grid of possible paths of freeze-out for S-bearing molecules. 


Model 

S—>-grains 

HS^grains 

HS 2 —>'grains 

S 2 ^grains 

CS-)-grains 

HCS^grains 

FI 

100% mS 

100 % mHS 

100 % mHS 2 

100 % mS 2 

100 % mCS 

100 % mHCS 

0 % mH 2 S 

0 % mHaS 

0 % mH 2 S 2 

0 % mHzS 

0% mHCS 

0 % mH 2 CS 

F2 

0% mS 

0% mHS 

0 % mHS 2 

0 % mS 2 

0% mCS 

0% mHCS 

100 % mH 2 S 

100 % mH 2 S 

100 % mH 2 S 2 

100 % mH 2 S 

100% mHCS 

100 % mH 2 CS 

F3 

50% mS 

50% mHS 

50% mHS 2 

50% mS 2 

50% mCS 

50% mHCS 

50% mH 2 S 

50% mHaS 

50% mH 2 S 2 

50% mH 2 S 

50% mHCS 

50% mH 2 CS 


sulfur residue (S-residue) that we are now able to produce 
on the grain. As we can see, in the updated version of the 
code, S-residue is in fact the fifth most abundant ice-mantle 
S-bearing species. These findings are extremely important, 
both because they match what has been alr e ady e xperimen- 
tally determined bv IJimenez-Escobar et aP (l2014h as an at¬ 
tempt to explain the lack of sulfur in dense regions. We now 
need to provide a better estimate of this residue in order to 
understand how its presence might affect the total amount 
of interstellar sulfur; we have decided not to speculate more 
about its nature, since only experiments can give us informa¬ 
tion on the chemical structure of these ref r actory molecules, 
but refer the interested reader to ISteudell (120031 1. 

Based on the computational work that has previously 
been performed in this area (see Introduction) as well as on 
some simple test models, we have not investigated changes 
in the temperature or in the density of our sources; instead, 
we analyse how the calculated abundances vary with the 
cosmic-ray ionisation rate by choosing a standard interstel¬ 
lar (Co = 1.3x 10 “^^ s“^), an enhanced (1.3x 10 “^® s“^) and 
a super-enhanced (1.3x10“^® s“^) value for this parameter. 
With these parameters we then repeat our sensitivity tests 
for each freeze-out chemistry (FI, F2, F3; Table[7|. Figure^ 
and Fig. |4] display the most abundant species found at the 
end of the collapse phase (where temperatures are as low as 
10 K and only non-thermal desorption can occur) and when 
the hot core is finally formed (^100 K), respectively. Note 
that while Fig. [3] shows the different molecules in the solid 
state (on the grain surface), in Fig. U we report the species 
in the gas phase, with the exception of the S-residue. Both 
figures show the phase where the bulk sulfur is located, and 
despite the increase in C, cosmic-ray induced desorption in 
not efficient enough in our models to return a significant 
amount of sulfur to the gas phase in Phase I. Specifically, 


in models with C = Co, we find about 0 . 02 % of sulfur in the 
gas phase. In models with C = lOOCo, this percentage rises 
to ~0.5% of the elemental sulfur in the gas phase at the end 
of Phase I. 

Phase I appears to be dominated by H 2 S (in agreement 
with the scheme of iHatchell et al1ll998l l. where the abun¬ 
dance seems to be influenced by the degree of ionisation as 
well as by the percentage of hydrogenation. Starting the dis¬ 
cussion of trends in the results by considering the case with 
0 per cent hydrogenation (top panels), going from left to 
the right in Figure [3] the increase in the C values (from the 
standard to the super-enhanced) leads to a more efficient 
non-thermal desorption of H 2 S, which therefore decreases in 
the icy mantle in favour of other S-bearing species. The im¬ 
mediate consequence is a greater amount of hydrogen sulfide 
in the gas phase. The latter species then reacts with atomic 
ions (C”*", S’*", H^) at early times and molecular ions (H 3 O’'', 
HCO’^) at late times to liberate atomic sulfur, which goes 
on to react with OH in order to produce SO and SO 2 . These 
oxides are then frozen back on to the grain. Moreover, we 
highlight a predictable increase in the S-residue abundance 
(since the rate of cosmic ray ionisation directly depends on 
the C, value) and an increasingly inefficient formation of solid 
OCS, which is mainly produced on the grain by reaction 
between H 2 S and CO. Moving to the bottom panels, the 
chemical behaviour is altered by the higher level of hydro¬ 
genation. At the beginning the situation is very similar to 
the one described above, but once gaseous HS and S are 
formed, they will both freeze back on to the grain as H 2 S. 
Therefore, assuming a super-enhanced cosmic ray ionisation 
rate, compared to the previous case, we now observe a larger 
amount of H 2 S, OCS and S-residue in the solid state. 

Phase II is dominated by SO 2 and the second most 
abundant species is SO. This result is not surprising. If we 
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Cosmic ray enhancement (^o) 

Figure 3. The most abundant S-bearing species in the icy mantle at the end of the collapse phase (Phase I). 0, 50 and 100 on the vertical 
axis refer to the percentage of hydrogenation chosen (see text); 1, 10, 100 indicate a standard, an enhanced and a super-enhanced cosmic 
ionisation rate, respectively. 


look at the chain of reactions mentioned above, we actu¬ 
ally find that gaseous H 2 S, S and HS will eventually lead to 
the formation of SO and SO 2 . The chemistries of these two 
species are strongly related to each other: SO forms SO 2 by 
reaction with O or OH; SO 2 produces SO when reacting with 
C. We notice immediately from Fig.|3]that the SO/SO 2 ratio 
is sensitive to the cosmic-ray ionisation held, increasing with 
held strength. Finally, as we can see in Fig. [4] an observable 
abundance of S-residue is still locked inside the grain, even 
after the ice mantle has been completely desorbed. 


4.2 Comparison of theoretical results with recent 
observations 

In order to verify the reliability of our calculations, we looked 
at the observed abundances of 21 species in hot cores and 
high-mass SFRs reported in the literature in the last hve 
years. The reader should bear in mind that when performing 
a comparison such as this, order of magnitude agreement be¬ 
tween observation and model results is excellent. Telescope 
beams often cover large regions, resulting in average column 
densities and abundances (averaged over both temperature 
and density variations), and since it is observationally very 
hard to determine the evolutionary stage of a hot core, it be¬ 
comes equally hard to identify a ‘t = 0’ to which to anchor 
the model results. Despite these caveats, a comparison be¬ 


tween observational and model results is an important check 
of the model. 

iHerpin et al.l (l2009l ) surveyed four hot cores using the 
IRAM 30 m and CSO telescopes, and detected li nes of CS, 
PCS, H 2 S, SO, SO 2 and their ^"*8 isotopologues. foin et al.l 
1 I 2 OIOI) used the SMA interferometer to observe G19.61-0.23, 
detecting 17 molecules, including SO, SO 2 , OCS and CS 
and various isotopologues. Of these four molecules, only the 
SO abundanc e was determined using more than one line. 
iNeufeld et al.l ll2012l) discovered the mercapto radical (HS) 
in the high-mass SFR, W49N, and al so provided a m e asure - 
ment of the H 2 S abundance there. IZernickel et al.l (l2012l) 
performed a sub-mm survey of the high-mass SFR NGC 
63341, detecting over 20 molecules in two hot cores. Sev¬ 
eral of these molecu les are sulfur-bear ing: OCS, CS, SO, 
H 2 S, SO 2 and H 2 CS. IXU fc Wand ll2013l) also used the SMA 
interferometer to observe a massive SFR, G20.080.14N, de- 
tecting 11 molecula r species, including SO and SO 2 . Finally, 
iGerner et al.l (l2014l) observed a large sample of 59 high-mass 
SFRs at various stages of evolution, including 11 hot molec¬ 
ular cores. Sixteen molecules were detected, of which four - 
^®CS, SO, C®^S, OCS - contain sulfur. Note that all these 
cores span a range of masses, distances and evolutionary 
stages. 

We compared fractional abundances calculated by our 
models at 10® yr of Phase H, for three values of the cosmic- 
ray ionisation rate in Fig. [5] Variations in hydrogenation 
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Figure 4. The most abundant S-bearing species in the gas (excluding the residue, which is solid) at the end of the desorption phase 
(Phase II). 0, 50 and 100 on the vertical axis refer to the percentage of hydrogenation chosen (see text and Table (Tj; 1, 10, 100 indicate 
a standard, an enhanced and a super-enhanced cosmic ray ionisation rate, respectively. 


efficiency during Phase I have little effect on the gas-phase 
abundances in Phase II, and thus we present results from 
“50 per cent hydrogenation” models only. In general, the 
agreement between calculated abundances and abundances 
derived from the observations is very good: models and ob¬ 
servational results are in agreement for two thirds of the 
species, and in the some of the remaining third, disagree¬ 
ment is less than an order of magnitude. There are issues 
with species which are formed mostly by UV and CR dis¬ 
sociation: CN, N 2 H+, C 2 H and HCO+ are under-estimated 
by our model. In the model, interstellar UV photons are as¬ 
sumed to become extinguished in the outer core material. 
Only cosmic-ray induced UV photons play an active part 
in the chemistry of our simulations. Hot cores, in reality, 
may be significantly fragmented, allowing UV photons to 
penetrate more deeply into the core than in a homogeneous 
medium. This may account for our under-estimation of these 
species in our model. There is little evidence in our galaxy 
for cosmic- ray ionisation r ates higher than those tested in 
our model (iDalgarndbOOd ). 

Observations indicate that CS is the most abundant 
sulfur-bearing species, however models struggle to produce 
enough CS to match observed amounts. This is a recur- 
rent problem in war m-temperature chemical models (e.g., 
IWakelam et ahlboill ). This under-estimation may stem from 
a lack of understanding of the warm-temperature chemistry 
of CS, incorrect rate coefficients for its gas-phase reactions. 


shock passage through the observed regions l|Viti et al.l 
l200ll) . etc. Other sulfur-bearing species are observed to have 
similar abundances to each other: OCS, H 2 S, SO, SO 2 , H 2 CS 
and HS, all have average fractional abundances in the range 
7-70x10“® in Fig. [5] Simulation results are within the same 
range for all the above species with the exception of SO 2 , 
which seems to be over-estimated, and SO, which appears 
to be under-produced (the relationship between these two 
species has already been discussed in detail - see Sect. 14.11) . 

In Table [21 we present the abundances of sulfur-bearing 
species for three values of C, the cosmic-ray ionisation rate. 
Note that we also include the abundances of species which 
have not been observed to date and the amount of a poten¬ 
tial residue. No ions are shown, since they are of low abun¬ 
dance (a;(ion)^10“^^), and do not contribute significantly 
to the sulfur budget. As previously mentioned, apart from 
CS and SO 2 , agreement is good between calculated abun¬ 
dances of sulfur-bearing species and their observed values. 
An interesting finding is that the abundance of H 2 S 2 reaches 
observable levels; therefore, the absence of its detection to 
date might be only due to the fact that this molecule has a 
small dipole moment (~1.2D). 

From Table [9] it is evident that the abundances of 
most sulfur-bearing species are influenced by C,, with some 
molecules only being significantly present in standard C, en¬ 
vironments (HS, NS, HCS). In order to make these changes 
more evident, we plot the fractional abundance of the S- 
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Table 9. Calculated fractional abundances of selected sulfur-bearing spe cies as a function of cosmic-ray ionisa t ion rate, compared t o 
their mean observed values towards variou s hot cores. Data are tak en from iHeroin et alJ l|200gh ; iQin et ^ l|2010l ') ; iNeufeld et al JSqII); 
IZernickel et alJ ll2012ll ; IXu fc Wand ll2Q13h and iGerner et alJ ll2014l^ . 


Species 

Standard Model 

Enhanced Model 

Super-enhanced Model 

Observed range 

SO 2 

1.3x10“'^ 

1.3x10-° 

7.2x10-'^ 

7.9x10-11-3.9x10-® 

so 

4.9x10-® 

1.5x10-® 

6.1x10-'^ 

8.5x10-11-1.9x10-'^ 

OCS 

2.9x10-® 

5.5x10-® 

2.6x10-1° 

2.0xl0-i°-4.0xl0-^ 

H 2 CS 

1.1x10-® 

3.9x10-® 

trace 

1.0x10-® 

H 2 S 

8.5x10-® 

trace 

1.0x10-1® 

7.3x10-11-8.3x10-® 

CS 

1.0x10-1° 

4.1x10-1® 

3.1x10-1° 

1.8xl0-i°-5.0xl0-^ 

HS 

1.2x10-® 

trace 

trace 

7.0x10-® 

NS 

1.2x10-® 

trace 

trace 

4.2x10-® 

HCS 

5.5x10-11 

trace 

trace 


S 

3.2x10-11 

1.9x10-11 

3.5x10-1° 


HS 2 

4.0x10-1° 

2.3x10-1° 

trace 


H 2 S 2 

1.7x10-® 

6.8x10-1° 

trace 


S 2 

2.4x10-® 

trace 

trace 


mS2 

trace 

1.5x10-1® 

7.9x10-1® 


mCS 2 

trace 

trace 

trace 


S-residue 

5.5x10-® 

2.3x10-® 

8.5x10-® 



Standard model 


Enhanced model 


Super-enhanced model 




3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 

log{Time) [yr] 


4.5 5.0 5.5 6.0 6.5 7.0 

log(Time) [yr] 


4.5 5.0 5.5 6.0 6.5 7.0 

log(Time) [yr] 


Figure 6. Results from the warm-up phase of the model, showing the effect of the cosmic-ray ionisation rate upon the fractional 
abundances of various species. The left panel is extended slightly to show the desorption events (those of t he (i) pure spec ies on the 
surface of the ice, (ii) monolayer on H 2 O ice, (Hi) volcano desorption, and (iv) co-desorption with H 2 O. See lViti et al.ll2004l . for a full 
description). The right panel shows species labels. This figure highlights the issues with comparing model results at a certain time to 
observational values, since fractional abundances can vary by several orders of magnitude. 


bearing species as a function of the time during the desorp¬ 
tion phase (Fig. [6)). In addition, we show the variation of 
H 2 O, HNC and CO, as standards for onr models. In gen¬ 
eral, abundant molecules, such as SO and OCS, are present 
at observational levels at all three strengths of and may be 
good tracers of (( if their observed abundances can be tightly 
constrained. SO 2 is also abundant at high levels throughout 
the tests, but does not vary much in abundance. 


A further consideration evident in Fig. is that the 
abundances of some species vary considerably with time, and 
as stated in the introduction to this section, this makes com¬ 
parisons with observations challenging. For instance, gas- 
phase H 2 S abundances drop in all scenarios by almost six 
orders of magnitude, but on different time-scales, and H 2 CS 
abundances by almost four orders of magnitude. Other 
species vary in abundance by an order of magnitude or so. 
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Observed vs. calculated abundances 
in hot cores and SFRs 



Figure 5. Observed (solid coloured lines) and ranges of cal¬ 
culated abundances from our models (hatched regions) for 
the s pecies reported in a number of rec e nt observational pa¬ 
pers jHerpin et alj 120091: iQin et alj I2OI0I: iNeufeld et alj I2OI2I : 
IZernickeTemTI^Ol^lXufc^^hnell^ni^ : lGerner‘^mi201^ . deal¬ 

ing with hot cores and high-mass SFRs. 

whereas species like SO 2 are present at their solid-phase 
abundances, i.e., their abundances are unaltered by gas- 
phase chemistry. 

If we now consider the refractory sulfur-bearing species, 
S 2 , CS 2 and the sulfur residue, it is evident that an appre¬ 
ciable amount (up to 6 per cent) of sulfur in our model 
is locked in form of this residue. It is hard to speculate 
about the chemical nature of this residue; recent experi¬ 
ments djimenez-Escobar et ahl |2014| ) seem to attribute its 
existence to the formation of sulfur polymers, such as Ss- 


There may, of course, be several ways to sink sulfur into 
this refractory residue, and we only consider a single route 
here, led by the experimental evidence. 

Finally, we are able to produce detectable abundances 
of species such as H 2 S and CS but only on the grain surface 
(Table [8]). This means that, while there is a good under¬ 
standing of the mantle sulfur chemistry, the rates of some 
gas-phase reactions need to be refined; in particular, the key 
step seems to be the high efficiency of SO and SO 2 formation 
from S and HS as reactants. 


5 CONCLUSIONS 

In summary, we have extended the ucl_CHEM chemical 
model to include new experimental and theoretical results on 
sulfur chemistry; in particular, we have inserted laboratory 
data from experiments where H 2 S ice is processed by pro- 
tons which simulate c osmic ray impacts on the grain surface. 
iGarozzo et al. ll l20iB detected in their icy mantle analogues 
the formation of CS 2 , which has not been observationally de¬ 
tected in the gas-phase ISM yet. In our modelling, we find 
that CS 2 was only present in small amou nts in icy man¬ 
tles. Follow i ng on from the experiment by IGarozzo et al.1 . 
IWard et al.1 (|2ni2l ) studied the production of solid OCS from 
GS 2 and O ice. We therefore also added the latter reaction 
to our chemical network. OGS abundan ces were enhan ced 
by a factor of ^2, driven mainly by the IGarozzo et al.1 re¬ 
action, mH 2 S -I- mGO —» mOCS - I- H 2 , and only partially 
by the reaction from IWard et ahl (l2012fl . Finally, we also 
revised our gas-pha se re actions based o n recen t papers by 
iLoison et al.r(l2012f) and lAdriaens et al.1 ll20ld) . We looked 
at the influence of the new reaction network on the fractional 
abundances of selected S-bearing species by comparing the 
results with the original version of our code. Abundances of 
species such as HS 2 and H 2 S 2 increased by large factors in 
the solid phase via very efficient hydrogenation of accreting 
species, but not significantly enough to make them major 
carriers of sulfur. Our main result was that a residue, left 
on the grain surface after ice desorption, could harbour a 
large a mount of sulfur (10~^), in l ine with recent experi¬ 
ments djimenez-Escobar et al. l2014h . 

A comparison with astronomical observations was also 
carried out. A good agreement is obtained although we find 
it difficult to reconcile theoretical estimates with observa¬ 
tions of some of the most abundant species such as H 2 S, 
GS and OGS. Our results show that at early stages during 
the collapse phase of star-formation, H 2 S is the predominant 
sulfur-bearing species in the icy mantles. As the evolution 
of a pre-stellar core continues SO 2 and SO are efficiently 
formed, in agreement with earlier work. In this respect, the 
pivotal processes appear to be the gas phase production of 
SO and SO 2 from S and HS as reactants. Another key fac¬ 
tor to be considered is the presence of a sulfur residue on 
the grain surface, which may affect the observability of some 
species in the gas-phase. This present paper is therefore an 
important step forward in understanding the sulfur chem¬ 
istry in regions of star formation as we constrain the problem 
to two main factors; the gaseous interactions between S or 
HS and O or OH to form SO and SO 2 and the existence of 
a sulfur residue on the grain surface. 
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